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Abstract 

We show that a gas thermometer in contact with a stationary classical system 
out of thermal (Boltzmann) equilibrium evolves, under very general conditions, to- 
wards a state characterized by a Levy velocity distribution. Our approach is based 
on a kinetic-like equation that applies to a wide class of models for the system- 
thermometer interaction. The results clarify the role of non-exponential energy dis- 
tributions as possible generalizations of the Boltzmann distribution for systems 
where the usual formulation of thermostatistics may not apply. In particular, they 
show that the power-law distributions derived from Tsallis's nonextensive formal- 
ism are irrelevant to the stationary state of the thermometer, thus failing to give 
a consistent description of the system-thermometer equilibrium. We point out the 
need of a generalized thermostatistical formulation able to give a unified frame to 
Levy and Maxwell distributions. 
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Levy distributions [1] have been extensively applied to the description and 
modeling of a wide class of physical processes, ranging from anomalous trans- 
port in disordered media [2] and turbulent flows [3], to phase-space diffusion 
in dynamical systems [4] and polymer dynamics [5]. They have also found 
application in other branches of science, such as in biology [6,7]. The ubiq- 
uity of Levy distributions in many natural phenomena is a straightforward 
consequence of their stable character under summation of random variables 
whose distributions have diverging moments, as shown by Paul Levy in his 
generahzation of the central hmit theorem [8] . This essential result gives Levy 
distributions the same status as the Gaussian distribution in the statistical 
description of stochastic processes. 

While most applications of Levy distributions deal with dynamical and trans- 
port processes, we show in this Letter that they replace the equilibrium Maxwell 
distribution — a Gaussian in the velocity variable — in an extended version of 
the scenario of a thermal measurement. Wc analyze the asymptotic energy 
distribution of a thermometer in thermal contact with a classical system in a 
nonequilibrium stationary state. While the nature of this nonequilibrium state 
is not explicitly specified, we typically refer to a situation where the system is 
not isolated and suffers the effect of external stationary forces that maintain its 
energy distribution apart from the prediction of Boltzmann-Gibbs statistics. 
It is an essential fact of thermostatistics [9] that, if the system is in thermo- 
dynamical equilibrium, its Boltzmann energy distribution is "copied" by the 
thermometer in such a way that its average energy per degree of freedom, i.e. 
its temperature, becomes equal to that of the system. Now, what properties 
of the state of the system are detected by the thermometer in the case that 
the system is maintained out of equilibrium? 

The motivation of this question is two- fold. First, since all physical systems 
are to some extent thermodynamically open and subject to external influence, 
the possible effect of these factors in a thermal measurement constitutes a 
problem of empirical relevance. In general, this problem transcends the lim- 
its of equilibrium thermostatistics and calls for a nonequilibrium, dynamical 
formulation. We choose a kinetic-like approach which describes a wide class 
of interaction models for the thermal contact between the system and a gas 
thermometer. Our results show that the thermometer attains a velocity distri- 
bution which is fully determined by the high-energy distribution of the system, 
and results to be insensitive to other details in the system state. In the case 
where the energy distribution of the system decays as a power law, as ad- 
vanced above, the thermometer velocities approach a Levy distribution. This 
distribution "copies" the singularities associated with divergent energy mo- 
ments in the system, and reduces to a Maxwellian when the average energy is 
flnite. 

The second motivation for our question has to do with the possibility of gen- 



2 



eralizing thermostatistics to the description of physical systems where the 
usual Boltzmann-Gibbs formulation may not apply [10]. Tsallis's nonextensive 
thermostatistics, for instance, is claimed to replace the usual formulation for 
systems where long-range interactions lead extensivity assumptions to break 
down, while containing Boltzmann-Gibbs statistics as a limiting case [11]. In 
Tsallis's formalism, the Boltzmann exponential distribution is replaced by a 
class of power-law functions derived from a variational principle for a gener- 
alized entropy. Whether this formalism is compatible with the basic facts of 
equilibrium between systems in thermal contact can be decided, precisely, by 
studying the interaction of a thermometer and a system with such power-law 
energy distribution. Our results show that, since Levy distributions are not 
regarded in this formalism as possible equilibrium states, Tsallis's thermo- 
statistics fails to give a consistent picture of thermo dynamical equilibrium. 

Consider a system in a stationary state which does not necessarily coincide 
with Boltzmann equilibrium. We assume that the state of the system is macro- 
scopically characterized by an energy distribution Fo{e). The system interacts 
with a thermometer consisting of an ensemble of independent particles of mass 
m moving in a d-dimensional domain. Within this domain, the distribution of 
particles is supposed to be spatially homogeneous and isotropic in velocities, 
such that the associated distribution function F{e, t) depends on the energy 
e and the time t only. The interaction between system and thermometer is 
assumed to fulfill the following conditions, (i) Interaction events are time lo- 
calized and infrequent, such that the typical time between interactions is large 
as compared with the relaxation time of the energy distribution of the sys- 
tem toward Fo(e). This insures that (a) at each event, the system is found 
at its stationary state and (b) any correlation between the states of system 
and thermometer, created by the interaction, dies out before the next event 
takes place, (ii) The amount of energy interchanged at each event is small as 
compared with the total energy of both the system and the thermometer. 

Within these conditions, it is possible to describe the evolution of F{e, t) by 
means of a kinetic-like equation. On the other hand, as a consequence of con- 
dition (i), no evolution equation is required for the system distribution. To 
represent the interaction between system and thermometer we choose a very 
generic class of models — closely related to stochastic Maxwell models [12] — 
which considerably simplifies the mathematical treatment of the kinetic prob- 
lem and, at the same time, does not imply severe limitations to the physics 
of the interaction. In fact, these models include interactions where momen- 
tum conservation can either hold or be violated, and admit arbitrary angular 
dependence in the associated cross sections. They assume that the interac- 
tion rate is independent of the energy and that energy itself is conserved. The 
corresponding evolution equation for F(e, t) stands for the balance between 
events where, due to interaction with the system, the thermometer reaches or 
abandons its states of energy e. Such events yield, respectively, positive and 
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negative contributions to dtF. The equation is more conveniently written in 
the Fourier representation, for the Fourier-transformed distribution 

= = [ exp{-ik-v)f{v,t)dv'^, (1) 



with u = k"^ /2m. Here, 



, T{d/2) 2 d 



F{mv^/2,t) 



(2) 



is the (isotropic) velocity distribution associated with F{e,t) [12]. Due to the 
normalization of -F(e), we have $(0,t) = 1 for all t. With an analogous defi- 
nition for ^o{u) in terms of Fo{e), the Fourier-transformed evolution equation 
reads [12]: 

1 

dt^u, t)^ j w{s)^q{su)^{1 - s)u, t\ds - ^{u, t). (3) 





The kernel w{s) characterizes the interaction, and satisfies the normalization 
Jq w{s)ds — 1. When necessary, we assume that the moments 

1 

= / s^w{s)ds (4) 



are well defined. 



The above condition (ii) imposes that w{s) is appreciably different from zero 
only for s 0. Taking this condition into account, we write ^[{1 — s)u,t] ^ 
^{u,t) — sudu^{u,t), to obtain an approximate form of Eq. (3). Its solution 
can be shown to asymptotically approach the stationary distribution 



^{u) — exp 



1 rdv' 

Wi J u' 




w{s) ln^o{su')ds 



(5) 



Since ~ except for s ~ the asymptotic solution is fully deter- 
mined by the behavior of $o(w) close to u — 0. We pay particular attention 
to the special case where, for u ^ 0, 

$o(m) ~ 1 - aou", (6) 



with ao > and < u < 1. For u ^ 1, this form corresponds to an energy 
distribution with a power- law tail for high energies, Fo{e) ~ e"*^""^. While Fq 



4 



is still normalized, its first moment — namely, the average energy — diverges. In 
the limit z/ = 1 a finite average energy is recovered and the long high-energy 
tail is lost. For ^o{u) as in Eq. (6), the stationary Fourier distribution for the 
thermometer is 



with a = uqW^I vwi. Note that, for < z/ < 1, reproduces the sin- 

gularity of $o(w) at M = 0. The corresponding Fourier-transformed velocity 
distribution, (j){k) = $(A;^/2m) [see Eq. (1)], reads 



with b — a/{2my and 7 = 2i/ (0 < 7 < 2). Note that, in the velocity 
representation, this corresponds to a Levy distribution for each velocity com- 
ponent. In fact, the Fourier transform of, say, f{vx) = J dvydv^ ■ ■ ■ f{v) is 
(j){kx, 0,0, ■ ■ ■) = exp(— 6|A;^|'^). Though, in general, it is not possible to write 
explicitly the corresponding form of f{vx), it is known that for large \vx\ and 
7 7^ 2 it behaves as f{vx) ~ l^xl"'''"^ [!]■ For 7 = 2, f{vx) is a Gaussian 
distribution. 

We stress that the result (8) for the Fourier-transformed velocity distribution 
of the thermometer is independent of the detailed form of the distribution 
Fo(e) of the system at low energies. As long as condition (ii) is satisfied, 4>{k) 
is completely determined by the power- law high-energy tail of Fo(e). 

Let us analyze these results in the few cases where the velocity distribution 
can be explicitly obtained. We first consider the situation where the average 
energy of the system under study is finite, E = J eF(){e)de < 00. Under this 
condition we have, in Eq. (6), = 2E/d and u = 1. This gives, for the 
Fourier-transformed velocity distribution of the thermometer, the Gaussian 
(f){k) — e:K.p{—Ek'^/md). In the velocity representation we get the Maxwellian 



The average value of the kinetic energy over this distribution is, as may be 
expected, {mv'^/2) = E. In other words, the stationary velocity distribution of 
the thermometer is a Maxwellian whose average energy per particle coincides 
with that of the system. This result includes of course the case where the 
system is itself in thermodynamical (Boltzmann) equilibrium at temperature 
T, for which E ocT. 

For u — 1/2, the energy distribution of the system decays as -Fo(e) ~ e~^/^ and. 



(7) 



=exp(-6F), 



(8) 




(9) 
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consequently, the average energy diverges. The Fourier-transformed velocity 
distribution of the thermometer results to be (^{k) = exp{—bk), with b = 



(y.QWii2\J'2,/mwi. Taking into account its units, we write b = po/m, where Pq 
is a characteristic linear momentum given by the energy distribution Fo{e) 
and the specific form of the interaction kernel w{s). In velocity space, each 
component Vi is distributed according to a Cauchy distribution, 

f{vi) = 2^ T^- 10 

TT Po + ITT' "^i 



Note that, here, the relevant dynamical variables in the thermometer dis- 
tribution are the momentum components mvi and, accordingly, the relevant 
parameter is the momentum po- Compare with the result for u = 1, Eq. (9), 
where the distribution is a function of the energy and the relevant parameter 
is E. 

The present results shed light on the role of non-Boltzmannian energy dis- 
tributions in defining the stationary states of out- of- equilibrium interacting 
systems. In particular, they show that, under the above conditions (i) and 
(ii), the Gaussian stationary velocity distribution of an ensemble of indepen- 
dent particles interacting with a non-equilibrium stationary system results to 
be a special instance of the more general situation where the velocity compo- 
nents have Levy distributions. We stress that, despite the fact this conclusion 
is limited by the validity of the assumptions (i) and (ii) , these conditions are 
fully compatible with the paradigm of thermal measurement [9] and therefore 
refer to a realistic, experimentally accessible situation. 

As advanced above, our conclusions are particularly significant in the evalua- 
tion of generalized formulations of thermostatistics as plausible descriptions of 
physical systems to which the usual Boltzmann-Gibbs theory may not apply. 
We focus our attention on Tsallis's formulation, which has motivated large 
amounts of work over the last two decades [10]. Within Tsallis's statistics, 
canonical probability distributions are obtained from maximization of a gen- 
eralized entropy [11] 

Sg-- ^~^f , (11) 



where pi is the probability of state i, and the sum runs over all the states 
of the system. The parameter q quantifies the nonadditivity of Sg, and is 
therefore a measure of nonextensivity in the system under consideration. The 
maximization process is subject to the generalized constraint [13] 

Eg = (12) 
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with ei the energy of state i, and leads to the probabihty distribution 

^p(60 = Z-'[l - (1 - q)(5,{e, - E,)f^'-'^\ (13) 

where /5q is a ( "renormahzed" [13]) Lagrange multiplier, analogous to the 
inverse temperature, and Zq is a normalization constant, analogous to the 
partition function. Boltzmann-Gibbs formulation is recovered for — > 1. Note 
that, for systems whose density of states behave as p(e) oc e'^ — a class which 
includes most classical systems — the probability p{t) corresponds to a Fourier- 
transformed energy distribution of the form of Eq. (6), with v = —a + {q — 
2)/(l-?). 

An attractive feature of Tsallis's generalized thermostatistics is that it pre- 
serves most of the mathematical structure of the usual theory, including the 
definition of thermodynamic functions, Legendre transformations, and even 
linear nonequilibrium properties [10]. In fact, a large part of the literature 
on this topic is devoted to the formal extension of thermodynamical relations 
to the generalized formulation. Thus, putting aside its phenomenological ap- 
plications to nonequilibrium processes such as anomalous diffusion [14] and 
turbulence [15], the real tour de force of Tsallis's formalism lies in the de- 
scription of thermal equilibrium for systems where Boltzmann-Gibbs theory 
is supposed to fail. This is in fact the original and most frequently invoked mo- 
tivation of the formulation and, consequently, constitutes its genuine source of 
validation. Our results clarify whether Tsallis's statistics is relevant to an as- 
pect that — in spite of its essential role in thermodynamics — has been scarcely 
treated in the profuse literature on Tsallis's theory [16], namely, thermal equi- 
librium between interacting systems. 

Assume to have a system which, due to its nonextensive nature, exhibits an 
energy distribution of the form of Eq. (13), as predicted by Tsallis's thermo- 
statistics. Allow furthermore the system to interact with a thermometer, as 
specified above. Consider first that the nonextensivity index q and the density 
of states of the system are such that E = J2i ^iPi < oo. According to our re- 
sults, the stationary velocity distribution of the thermometer is a Maxwellian, 
corresponding to a Boltzmann exponential probability for the energy. As ex- 
pected for an ensemble of independent particles, the thermometer behaves 
as an extensive system {q = 1). The value of the parameters in the ther- 
mometer distribution predicted by Tsallis's formalism, however, results to be 
wrong. Abe and Rajagopal [16] have shown that a formal extension of equilib- 
rium conditions for two systems in thermal contact yields, in Tsallis's theory, 
Pq — Pqi, where q and q' are the nonextensivity indices corresponding to the 
two systems [17]. Our results show that this relation is generally not satisfied. 
In fact, taking g' = 1 as the nonextensivity index of the thermometer, Eq. (9) 
implies (3q' oc E~^. On the other hand, (3q shows in general a more complicated 
dependence on E. For instance, it can be readily shown that, for a system with 
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a density of states p(e) oc e'^, one has Pq oc Let us also point out 

that the generahzed average Eq of Eq. (12) plays no role in the parameter that 
defines the equilibrium distribution of the thermometer — namely, the inverse 
temperature (3qi — in spite of the fact that this generalized average replaces E 
in the formal extension of thermodynamics to Tsallis's formalism. 

In the case that the average Yl,i ^iPi diverges, our results imply that the ther- 
mometer approaches an energy distribution not belonging to the class of the 
system distribution, Eq. (13). Thus, no well-defined index q can be assigned 
to the thermometer. In this situation, the assumption that Tsallis's thermo- 
statistics would describe the system-thermometer equilibrium fails drastically. 
The only trace of the energy distribution of the system in the thermometer 
distribution arises from its power-law tail, which defines the Levy exponent 
7 = 2i/ in Eq. (8). Note that, again, the generahzed average Eq is irrelevant 
to the determination of the equilibrium distribution of the thermometer. 

Summing up, our kinetic description of a thermometer in contact with a sys- 
tem in a stationary state different from Boltzmann equilibrium, suggests that 
a suitable extension of thermostatistics should yield Levy distributions as the 
generalization of the Maxwellian velocity distribution, instead of the power- 
like functions of Eq. (13). We conclude that the energy probabilities derived 
from Tsallis's thermostatistics do not play the role of equilibrium distributions 
for systems in thermal contact. Though this conclusion is not completely gen- 
eral, it applies to a wide class of interaction models in the realistic situation 
where a system is put in contact with a thermometer. The qiicstion remains 
open as to which form of the entropy and which constraints should be used to 
derive Levy distributions from a variational principle. 

We thank G. Abramson and L Samengo for their critical reading of the 
manuscript. 
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